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Abstract 

It is shown, that the BOLTZMANN-like equations allow the formulation of a very 
general model for behavioral changes. This model takes into account spontaneous 
(or externally induced) behavioral changes and behavioral changes by pair interac- 
tions. As most important social pair interactions imitative and avoidance processes 
are distinguished. The resulting model turns out to include as special cases many 
theoretical concepts of the social sciences. 

A Kramers-Moyal expansion of the BOLTZMANN-like equations leads to the 
Boltzmann-Fokker-Planck equations, which allows the introduction of "social 
forces" and "social fields" . A social field reflects the influence of the public opinion, 
social norms and trends on behaviorial changes. It is not only given by external 
factors (the environment) but also by the interactions of the individuals. Variations 
of the individual behavior are taken into account by diffusion coefficients. 
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1 Introduction 



The methods of statistical physics have shown to be very fruitful in physics, but in the 
last decades they also have become increasingly important in interdisciplinary research. 
For example, the master equation has found many applications in thermodynamics 0, 
chemical kinetics [fj], laser theory |J and biology Moreover, Weidlich and Haag 
have successfully introduced it for the description of social processes || || like opinion 
formation M, migration [Q, agglomeration M and settlement processes [|TD| . 

Another kind of wide-spread equations are the Boltzmann equations, which have been 
developed for the description of the kinetics of gases [TT| and of chemical reactions |T2[ . 



However, BoLTZMANN-/zA;e equations |L3| play also an important role for quantitative 
models in the social sciences: It turns out (cf. sect. ^2|) that the logistic equation for 
the description of limited growth processes [|14|, [L5|], the socalled gravity model for spatial 



exchange processes [16|, and the game dynamical equations modelling competition and 



cooperation processes [|H], [18| are special cases of BoLTZMANN-like equations. Moreover, 
BoLTZMANN-like models have recently been suggested for avoidance processes of pedestri- 
T9L EU[ and for attitude formation by direct pair interactions of individuals occuring 



ans 



in discussions [19, 21 



In this paper we shall show that BoLTZMANN-like equations and Boltzmann-Fokker- 
Planck equations |H| are suited as a foundation of quantitative behavioral models. For 



this purpose, we shall proceed in the following way: In section ^| the BoLTZMANN-like 
equations will be introduced and applied to the description of behavioral changes. The 
model includes spontaneous (or externally induced) behavioral changes and behavioral 
changes by pair interactions of individuals. These changes are described by transition 
rates. They reflect the results of mental and psychical processes, which could be simulated 



with help of Osgood and Tannenbaum's congruity principle [E2|, Heider's balance 



theory |23| or Festinger's dissonance theory ||24| . However, it is sufficient for our model 



to determine the transition rates empirically (sect. |5|). The ansatz used for the transition 
rates distinguishes imitative and avoidance processes, and assumes utility maximization 
of the individuals (sect. |2.1|) . It is shown, that the resulting BoLTZMANN-like model for 
imitative processes implies as special cases many generally accepted theoretical approaches 
in the social sciences (sect. j2~2l ). 

In section ^ a consequent mathematical formulation related to an idea of Lewin 
is developed, according to which the behavior of individuals is guided by a social field. 
This formulation is achieved by a Kramers-Moyal expansion of the BoLTZMANN-like 
equations leading to a kind of diffusion equations: the socalled Boltzmann-Fokker- 



Planck equations [13]. In these equations the most probable behavioral change is given 
by a vectorial quantity that can be interpreted as social force (sect. |3.1| ). The social force 
results from external influences (the environment) as well as from individual interactions. 
In special cases the social force is the gradient of a potential. This potential reflects the 
public opinion, social norms and trends, and will be called the social field. By diffusion 
coefficients an individual variation of the behavior (the "freedom of will") is taken into 
account. In section ^ representative cases are illustrated by computer simulations. 

The Boltzmann-Fokker-Planck modell for the behavior of individuals under the in- 
fluence of a social field shows some analogies with the physical model for the behavior of 
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electrons in an electric field (e.g. of an atomic nucleus) |13|, [19| (cf. Hartree's selfconsis- 
tent field ansatz Especially, individuals and electrons influence the concrete form of 
the effective social resp. electric field. However, the behavior of electrons is governed by a 
different equation: the SCHRODINGER equation. 

In physics, the Boltzmann-Fokker-Planck equations can be used for the description 
of diffusion processes p? . 



2 The BOLTZMANN-like equations 

Let us consider a system consisting of a great number N ^> 1 of subsystems. These 
subsystems are in one state x of several possible states combined in the set Q. 

Due to fluctuations one cannot expect a deterministic theory for the temporal change 
dx/dt of the state x(t) to be realistic. However, one can construct a stochastic model for 
the change of the probability distribution P(x,t) of states x(t) within the given system 
(P(x,t) > 0, ^ P(x,t) = 1). By introducing an index a we may distinguish A different 

A 

types a of subsystems. If N a denotes the number of subsystems of type a, we have ^ N a = 

o=l 

N, and the following relation holds: 

A N 

P(x,t) = J2^Pa(x,t). (1) 

a=l iV 

Our goal is now to find a suitable equation for the probability distribution P a (x,t) of 
states for subsystems of type a (P a (x,t) > 0, 2J Pa(x,t) = 1). If we neglect memory 



effects (cf. sect. |6.1|) , the desired equation has the form of a master equation |13, 19| : 



^-P a (x,t)= [w a (x\x';t)P a (x',t)-w a (x'\x;t)P a (x,t)]. (2) 



dt 

x' £S! 
(x'j^x) 

w a (x'\x;t) is the effective transition rate from state x to x' and takes into account the 
fluctuations. Restricting the model to spontaneous (or externally induced) transitions and 
transitions due to pair interactions, we have [|T^, [HJ: 

A 

w a (x'\x;t) :=w a (x'\x;t)+Y,Y, E N b w ab {x',y'\x,y;t)P b {y,t). (3) 

b=ly£tty'en 

w a (x'\x; t) describes the rate of spontaneous (resp. externally induced) transitions from x 
to x' for subsystems of type a. w a b(x', y'\x, y; t) is the transition rate for two subsystems 
of types a and b to change their states from x and y to x' and y' due to pair interactions. 



Inserting @ into (0), we now obtain the socalled BoLTZMANN-/i£;e equations [13, [19 

^-P a (x,t) = X) \w a {x\x';t)P a {x',t) -w a {x'\x;t)P a {x,t)] (4a) 
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A 

+ E E E E w ab (x,y'\x',y;t)P b (y,t)P a (x',t) 

6=1 x' en yen y> en 

A 

- E E E E w °b(x', y'\*, y; tWv, t)p a (x, t) (4b) 

b=i x' en yen y ' en 

with 

w ab (x',y'\x,y;t) := N b w ab (x' ,y'\x,y;t) . (5) 

Obviously, (Hb) depends nonlinearly on the probability distributions P a (x,t), which is 
due to the interaction processes. 

Neglecting spontaneous transitions (i.e., w a (x'\x;t) = 0) the BoLTZMANN-like equations 
agree with the Boltzmann equations, that originally have been developed for the de- 
scription of the kinetics of gases |TT[. A more detailed discussion can be found in ]TP[ . 



In order to apply the BOLTZMANN-like equations to behavioral changes we have now to 
take the following specifications given by table 1 (cf. flI5| ): 



STATISTICAL PHYSICS 


BEHAVIORAL MODELS 


system 


population 


subsystems 


individuals 


states 


behaviors (concerning a 
special topic of interest) 


types of subsystems 


types of behavior 
(subpopulations) 


transitions 


behavioral changes 


interactions 


imitative processes, 
avoidance processes 


fluctuations 


"freedom of will" 



Table 1: Specification of the notions used in statistical physics for an application to be- 
havioral models. 

It is possible to generalize the resulting behavioral model to simultaneous interactions of 
an arbitrary number of individuals (i.e., higher order interactions) [13j |I9| . However, in 
most cases behavioral changes are dominated by pair interactions. Many of the phenomena 
occuring in social interaction processes can already be understood by the discussion of 
pair interactions. 

2.1 The form of the transition rates 

In the following we have to find a concrete form of the effective transition rates 

A 

w a (x'\x; t) = ^(aj'laj; t) + E E E W ab (x'\x, y; t)P b (y, t) , (6a) 

6=i yen y'en 
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Wab(x'\x,y;t) := ^ w ab (x' ,y'\x,y;t) (6b) 

(cf. (^D) that is suitable for the description of behavioral changes. W ab (x'\x,y;t) is the 
rate of pair interactions 

x' ^- x, (7) 

where an individual of subpopulation a changes the behavior from x to x' under the 
influence of an individual of subpopulation b showing the behavior y. There are only two 
important kinds of social pair interactions: 

x' < — x (x' x) (8a) 
x' ^- x (x'^x). (8b) 

Obviously, the interpretation of the above kinds of pair interactions is the following: 



• The interactions (|8]a) describe imitative processes (processes of persuasion), that 
means, the tendency to take over the behavior x' of another individual. 

• The interactions (§|b) describe avoidance processes, where an individual changes the 
behavior when meeting another individual showing the same behavior x. Processes 
of this kind are known as aversive behavior, defiant behavior or snob effect. 



The corresponding transition rates are of the general form 

W ab (x'\x,y;t) := u ab (t)Rl b (x'\x; t)8 yx , (9a) 
+ u ab {t)R 2 ab (x'\x;t)5 yx , (9b) 

where the term ([|a) describes imitative processes and the term (Qi) describes avoidance 
processes. 5 yx has the meaning of the Kronecker function. By inserting (^) into (|5[) we 
arrive at the following general form of the effective transition rates: 

A 

w a (x'\x; t) := u a (t)R a (x'\x; t) + £ u ab (t) [Rl b (x'\x; t)P b (x', t) + R 2 ab (x'\x; t)P b (x, t) 

b=l 

(10a) 

For behavioral models one often assumes 

R k ab (x'\x;t) := f k ab (t)R a (x'\x;t) . (10b) 

In pi), 



• v a {t) is a measure for the rate of spontaneous (or externally induced) behavioral 
changes within subpopulation a. 

• R a (x'\x;t) [resp. R a (x'\x; t)} is the readiness for an individual of subpopulation a 
to change the behavior from x to x' spontaneously [resp. in pair interactions]. 

• v ab (t) = N b V ab (t) is the interaction rate of an individual of subpopulation a with 
individuals of subpopulation b. 
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• fab{t) * s a measure for the frequency of imitative processes. 

• fab{t) is a measure for the frequency of avoidance processes. 

A more detailled discussion of the different kinds of interaction processes and of ansatz 
(0) is given in||||. 



For R a (x'\x; t) we take the quite general form 

JJ a {x',t)-U a (x,t) 

R a (x'\x;t) = — — — (11a) 

1 1 ; D a (x',x;t) 1 ; 

with 

D a (x',x;t) = D a (x,x';t) > 

(cf. H |19|). Then, the readiness R a (x'\x; t) for an individual of subpopulation a to change 
the behavior from x to x' will be the greater, 

• the greater the difference of the utilities U a (.,t) of behaviors x' and x is, 

• the smaller the incompatibility ("distance") D a (x',x;t) between the behaviors x 
and x' is. 



Similar to Qllaj) we use 

e u a (x' ,t)-U a (x,t) 

"•fc>"> = m^j) ■ (llb) 

and, therefore, allow the utility function U a (x,t) for spontaneous (or externally induced) 
behavioral changes to differ from the utility function U a (x,t) for behavioral changes in 
pair interactions. Ansatz ( pT|) is related to the multinomial logit model |29|, |30| , and assumes 
utility maximization with incomplete information about the exact utility of a behavioral 
change from x to x', which is, therefore, estimated and stochastically varying (cf. |T9fl). 



Computer simulations of the BoLTZMANN-like equations (fj), (|10D, ([TT|) are discussed and 



illustrated in 1T9L ETL E3] (cf. also sect. 



2.2 Special fields of application in the social sciences 

The BoLTZMANN-like equations (0), ( |TUD include a variety of special cases, which have 
become very important in the social sciences: 



The logistic equation flbj, describes limited growth processes. Let us consider 
the situation of two behaviors x 6 {1,2} (i.e., P a (l,t) = 1 — P a (2,t)) and one 
subpopulation (A = 1). x = 2 may, for example, have the meaning to apply a 
certain strategy, and x = 1 not to do so. If only imitative processes 

2 ^- 1 (12) 
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and processes of spontaneous replacement 

1 < — 2 

are considered, one arrives at the logistic equation 



(13) 



jPxM = -z/ 1 (t) J R 1 (l|2;t)P 1 (2,t) + z/ 11 (t)/ 1 1 1 (t)i? 1 (2|l;t)(l-P 1 (2,t))P 1 (2,t) 

(14) 



= A(t)P 1 (2,t)(B(t)-P 1 (2,t 



The gravity model [0 describes processes of exchange between different places x. 
It results for R a (x'\x; t) = 0, f^(t) = 1, fl b {t) = 0, and A = 1: 



d_ 

It 



P{x,t) = v{t) £ 



JJ(x,t)-U(x',t) „U(x',t)-U(x,t) 



D(x,x' 



D(x',x) 



P(x,t)P(x',t) . 



(15) 



Here, we have dropped the index a because of a = 1. P(x,t) is the probability of 
being at place x. The absolute rate of exchange from x to x' is proportional to the 
probabilities P(x, t) and P(x', t) at the places x and x' . D{x, x') is often chosen as a 
function of the metric distance ||cc — cc'H between x and x'\ D(x,x') = D(\\x — x'\\). 

The behavioral model of Weidlich and Haag [§, ^|, || assumes spontaneous tran- 
sitions due to indirect interactions, which are, for example, induced by the media 
(TV, radio, or newspapers). We obtain this model for = = fl b {t) and 



U a (x, t) := 5 a (x, t) + J2 K ab Pb(x, t) . 

6=1 



(16) 



S a (x, t) is the preference of subpopulation a for behavior coupling param- 

eters describing the influence of the behaviorial distribution within subpopulation b 
on the behavior of subpoplation a. For > 0, K a b reflects the social pressure of 
behavioral majorities. 

The game dynamical equations [|I7|, [18], [19], result for /* 6 (t) = 5 a b, /^(i) = 0, and 



R a (x'\x; t) : = max (E a (x', t) - E a (x, t), 

(cf. |19], ^8|). Their explicit form is 
d 



dt 



P a (x,t) = [wa(x\x';t)P a (x',t) - w a (x'\x;t)P a (x,t) 



x'en 



+ v aa {t)P a {x,t) E a {x,t) - (E a ) 



(17) 

(18a) 
(18b) 



Whereas (|T^a) again describes spontaneous behavioral changes ( "mutations", inno- 
vations), fll8[b) reflects competition processes leading to a "selection" of behaviors 
with a success E a (x, t) that exceeds the average success 



(E a ) := ]T E a (x',t)P a (x',t). 

x'en 



(19) 
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The success E a (x,t) is connected with the socalled payoff matrices A ab = 



(Aabix, y)) by 



E a (x, t) := A a (x) + E E A*(aJ, y)P b (v, t) 

6=1 j/eC 



(20) 



rpTSl , |28|| . A a (a;) means the success of behavior a; with respect to the environment. 

Since the game dynamical equations fllBD agree with the selection mutation equations 
[17] they are not only a powerful tool in social sciences and economy |1| |32], |33| . 
but also in evolutionary biology [B3, 



3 The Boltzmann-Fokker-Planck equations 



We shall now assume the set Q of possible behaviors to build a continuous space. The n 
dimensions of this space correspond to different characteristic aspects of the considered 
behaviors. In the continuous formulation, the sums in (0), (^) have to be replaced by 
integrals: 



d_ 
~dt 



Pa(x,t) 



al n x' 



al n x' 



w a (x\x'; t)P a (x', t) - w a (x'\x; t)P a (x, t) 

w a [x'\x - x'\ t]P a (x - x', t) - w a [x'\x; t)P a (x, t)} , (21a) 



where 



w a [x'-x\x;t] :=w a (x'\x;t) := w a (x'\x; t) +J2 J <Pv J d n y f N b w ab (x',y'\x,y;t)P b (y,t) . 



(21b) 

A reformulation of the Boltzm ANN- like equations ( pTj) via a Kramers-Moyal expan- 
sion |38], |39|] (second order Taylor approximation) leads to a kind of diffusion equations: 



the socalled Boltzmann-Fokker-Planck equations 







d 



9 9 



Q aii (:z,t)P a (cM) (22a) 



with the effective dn/t coefficients 

K ai (x, t) := I oTx' {x\ - x l )w a (x'\x; t) 



(22b) 



and the effective diffusion coefficients^ 



Q aij (x,t) := Jd n x' (x'i -x i )(x' j - Xj)w a (x'\x;t) 
n 



(22c) 



Hn [ fL3| the expression for Q a ij(x,t) contains additional terms due to another derivation of (| 
However, they make no contributions, since they result in vanishing surface integrals (cf. p9[). 
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Whereas the drift coefficients K a i(x,t) govern the systematic change of the distribution 
P a (x, t), the diffusion coefficients Q a ij(x, t) describe the spread of the distribution P a (x, t) 
due to fluctuations resulting from the individual variation of behavioral changes. 

For ansatz (|10D , the effective drift and diffusion coefficients can be splitted into contribu- 
tions due to spontaneous (or externally induced) transitions (k = 0), imitative processes 
(k — 1), and avoidance processes (k — 2): 

K ai (x, t) = j2 K k ai (x, t) , Q aij (x, t) = £ Q k aij (x, t) , (23a) 

fc=0 fc=0 



where 



K° ai (x,t) := UaWjdTx'ix'i-xdRaix'lxit), 
K l Jx,t) := ^v^f^J^x'ix'i-xjR^x'lx^Ptix^t), 



6=1 

.4 



and 



K 2 ai (x,t) := J2"ab(t)f 2 ab (t) fcTx'ix'i-xjR^x'lx^P^t) (23b) 

6=1 J 

Q° aij (x,t) := u a (t)J<rxr(x? i -x i )tf j -x j )R a (x'\x-,t), 

Qli^t) ■= X>«*(*)/i(*) f^(x' i -x i )(a/ j -x j )R i {x , \x;t)P b {x , ,t) > 

6=1 J 

Q 2 aij (x,t) := SM*)/^*)/^^-^)^-^)^^;*)^*)- ( 23c ) 



The behavioral changes induced by the environment are included in K^ { (x,t) and 



3.1 Social force and social field 



The Boltzmann-Fokker-Planck equations 
equations (Langevin equations) 



are equivalent to the stochastic 



dt 



F ai (x,t) + Y,G mj (x,t)Z j (t) 



(24a) 



with 



and 



1 " 

K ai (x,t) = F ai (x,t) + - 

j,k=l 



_d_ 

dx k 



G a ij (x, t) 



Qui] ^ ' G qiki^C , t)(j a fcj(35, t) 
fc=l 



(24b) 



(24c) 
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(cf. [fL9|). For an individual of subpopulation a the vector £ a (x,t) with the components 



ta i (x,t) = 'EG aij (x,t)Zj(t) (25) 

describes the contribution to the change of behavior x that is caused by behavioral fluc- 
tuations £(t) (which are assumed to be delta-correlated and GAUSSian jl9"|). Since the 



diffusion coefficients Q a ij(x,t) and the coefficients G a ij{x } t) are usually small quantities, 
we have F a i(x,t) « K a i(x,t) (cf. (|24b ) ) , and ( 24a|) can be put into the form 



dx 

— ~ K a (x,t) + fluctuations. (26) 

Whereas the fluctuation term describes individual behavioral variations, the vectorial 
quantity K a (x,t) drives the systematic change of the behavior x(t) of individuals of 
subpopulation a. Therefore, it is justified to denote K a (x,t) as social force acting on 
individuals of subpopulation a. 

The social force influences the behavior of the individuals, but, conversely, due to inter- 
actions, the behavior of the individuals also influences the social force via the behavioral 
distributions P a (x,t) (cf. (pl|b), ( |22b| )). That means, K a (x,t) is a function of the social 
processes within the given population. 

Under the integrability conditions 

d d 
— K ai (x, t) = —K aj (x, t) for all i, j (27) 

OXj OXi 

there exists a time- dependent potential 

X 

V a (x,t):=-Jdx'-K a (x',t), (28) 

so that the social force is given by its gradient V: 

K a (x,t) = -WV a (x,t). (29) 

The potential V a (x,t) can be understood as social field. It reflects the social influences 
and interactions relevant for behavioral changes: the public opinion, trends, social norms, 
etc. 



3.2 Discussion of the concept of force 

Clearly, the social force is no force obeying the NEWTONian laws of mechanics. Instead, 
the social force K a (x,t) is a vectorial quantity with the following properties: 

• K a (x,t) drives the temporal change dx/dt of another vectorial quantity: the be- 
havior x(t) of an individual of subpopulation a. 
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• The component 

K ab (x,t) :=u ab (t) J d n x' (x' - x)[fl b (t)P b (x',t) + f 2 ab (t)P b (x,t)]R a (x'\x-,t) (30) 
n 

of the social force K a (x, t) describes the reaction of subpopulation a on the behav- 
ioral distribution within subpopulation b and usually differs from K ba (x,t), which 
describes the influence of subpopulation a on subpopulation b. 

• Neglecting fluctuations, the behavior x(t) does not change if K a (x,t) vanishes. 
K a (x,t) = corresponds to an extremum of the social field V a (x,t), because it 
means 

VV a (x,t) = 0. (31) 
We can now formulate our results in the following form related to Lewin's "field theory" 



25 : 



• Let us assume that an individuals objective is to behave in an optimal way with 
respect to the social field, that means, he or she tends to a behavior corresponding 
to a minimum of the social field. 

• If the behavior x does not agree with a minimum of the social field V a (x,t) this 
evokes a psychical tension (force ) 

K a (x,t) = -WV a (x,t) (32) 

that is given by the gradient of the social field V a (x,t). 

• The psychical tension K a (x,t) is a vectorial quantity that induces a behavioral 
change according to 

flip 

-^^K a (x,t). (33) 

• The behavioral change dx/dt drives the behavior x(t) towards a minimum x* a of the 
social field V a (x,t). 

• When the behavior has reached a minimum x* of the social field V a (x,t), it holds 

VV a (x,t) = (34) 
and, therefore, K a (x,t) = 0, that means, the psychical tension vanishes. 

• If the psychical tension K a (x,t) vanishes, except for fluctuations no behavioral 
changes take place — in accordance with (|33|) . 

In the special case, where an individual's objective is the behavior cc*, one would expect 
behavioral changes according to 

^ « 7 « - x) , (35) 
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which corresponds to a social field 



V a (x,t)^^(x* a -x) 2 (36) 

with a minimum at x*. Examples for this case are discussed in |40| . 

Note, that the social fields V a (x,t) of different subpopulations a usually have different 
minima x*. That means, individuals of different subpopulations a will normally feel dif- 
ferent psychical tensions K a (x,t). This shows the psychical tension K a (x,t) to be a 
"subjective" quantity. 



4 Computer simulations 

In the following, the Boltzmann-Fokker-Planck equations for behavioral changes will 
be illustrated by representative computer simulations. We shall examine the case of A = 2 
subpopulations, and situations for which the interesting aspect of the individual behavior 
can be described by a certain position x 6 [1/20, 1] on a one-dimensional continuous scale 
(i.e., n = 1, x = x). Then, the integrability conditions (|27|) are automatically fulfilled, 
and the social field 

X 

V a (x, t) = - f dx' K a (x', t) - c a (t) (37) 



is well-defined. The parameter c a (t) can be chosen arbitrarily. We will take for c a (t) the 
value that shifts the absolute minimum of V a (x, t) to zero, that means, 



/ X 

c a {t) := min I — / dx K a (x', 



t) . (3* 



• Since we will restrict the simulations to the case of imitative or avoidance processes, 
the shape of the social field V a (x, t) changes with time only due to changes of the 
probability distributions P a (x,t) (cf. (|23"1)), that means, due to behavioral changes 
of the individuals (see figures [1] to 

In the one-dimensional case one can find the formal stationary solution 

p -^ p -w||'( ! / tf w)' (39) 

\ x o / 

which we expect to be approached in the limit of large times t — > oo. Due to the depen- 
dence of K a (x) and Q a (x) on P a (x), equations (|39|) are only implicit equations. However, 
from (|39|) we can derive the following conclusions: 



If the diffusion coefficients are constant (Q a (x) = Q a ), Q3"9] ) simplifies to 

2 

"Qa 



P a (x) = P a (x ) exp (-^- \V a (x) + cj ) , (40) 



11 



that means, the stationary solution P a (x) is completely determined by the social 
field V a (x). Especially, P a (x) has its maxima at the positions x*, where the social 
field V a (x) has its minima (see fig. [1|). The diffusion constant Q a regulates the width 
of the behavioral distribution P a (x). 

• If the diffusion coefficients Q a (x) are varying functions of the position x, the behav- 
ioral distribution P a (x) is also influenced by the concrete form of Q a (x). From ( ]39| ) 
one expects high behavioral probabilities P a {x) where the diffusion coefficients Q a (x) 
are small (see fig. |2|, where the probability distribution P±(x) cannot be explained 
solely by the social field Vi(x)). 

• Since the stationary solution P a (x) depends on both, K a (x) and Q a (x), different 
combinations of K a (x) and Q a (x) can lead to the same probability distribution 
P a {x) (see fig. |] in the limit of large times). 

For the following simulations, we shall use the ansatz 

e U a (x',t)-U a (x,t) 

R a (x'\x; t) = — (41a) 

V 1 ; D a (x',x;t) V ' 

for the readiness R a (x'\x; t) to change from x to x' (cf. ( |lla|) ). With the utility function 

U °^-—K^ir)'- l - L fo < 41b > 

subpopulation a prefers behavior x a . L a means the indifference of subpopulation a with 
respect to variations of the position x. Moreover, we take 

Uab ® ■= e -W-*\/r r = E ui c ) 
D a (x',x;t) ' 20 ' K J 

where R can be interpreted as measure for the range of interaction. According to ( |41c|) . 
the rate of behavioral changes is the smaller the greater they are. Only small changes of 
the position (i.e., between neighboring positions) contribute with an appreciable rate. 



4.1 Sympathy and interaction frequency 

Let s a b(t) be the degree of sympathy which individuals of subpopulation a feel towards 
individuals of subpopulation b. Then, one expects the following: Whereas the frequency 
fab(t) of imitative processes will be increasing with s a b(t), the frequency fl b {t) of avoidance 
processes will be decreasing with s a b(t). This functional relationship can, for example, be 
described by 

fab(t) ■= ft(t) 8ab (t), 

fl{t) ■= fa(t)(l-s ab (t)) (42) 

with 

< s ab (t) < 1 . (43) 

fa(t) is a measure for the frequency of imitative processes within subpopulation a, flit) 
a measure for the frequency of avoidance processes. If we assume the sympathy between 
individuals of the same subpopulation to be be maximal, we have su(i) = 1 = S22(t). 
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4.2 Imitative processes (/*(£) = 1, flit) = 0) 



In the following simulations of imitative processes we assume the prefered positions to be 
xi = 6/20 and x 2 = 15/20. With 

(««*(*)) = (/iW) == ( I \ ) , (44) 

the individuals of subpopulation a = 1 like the individuals of subpopulation a = 2, but 
not the other way round. That means, subpopulation 2 influences subpopulation 1, but 
not vice versa. One could say, the individuals of subpopulation 2 act as trendsetters. 

As expected, in both behavioral distributions P a (x,t) there appears a maximum around 
the prefered behavior x a . In addition, due to imitative processes of subpopulation 1, a 
second maximum of Pi(x, t) develops around the prefered behavior x 2 of the trendsetters. 
This second maximum is small, if the indifference L\ of subpopulation 1 with respect to 
variations of the position x is low (see fig. Hj). For high values of the indifference L\ even 
the majority of individuals of subpopulation 1 imitates the behavior of the trendsetters 
(see fig. D! 

We shall now consider the case 

(*«*(*)) = (&(*)) == ( J I ) , (45) 

for which the subpopulations influence each other mutually with equal strengths. If the 
indifference L a with respect to changes of the position x is small in both subpopulations a, 
each probability distribution P a (x,t) has two maxima. The higher maximum is located 
around the prefered position x a . A second maximum can be found around the position 
prefered in the other subpopulation. It is the higher, the greater the indifference L a is (see 
fig. |). 

However, if L a exceeds a certain value in at least one subpopulation, only one maximum 
develops in each behavioral distribution P a (x, t)\ Despite the fact, that the social fields 
V a (x, t) and diffusion coefficients Q a (x,t) of the subpopulations a are different because of 
their different prefered positions x a (and different utility functions U a (x, t)), the behavioral 
distributions P a (x, t) agree after some time! Especially, the maxima x* a of the distributions 
P a (x,t) are located at the same position x* in both subpopulations. x* is nearer to the 
position x a of the subpopulation a with the lower indifference L a (see fig. |). 

4.3 Avoidance processes = 0, fl(t) = 1) 

For the simulation of avoidance processes we assume with x\ = 9/20 and x 2 = 12/20 that 
both subpopulations nearly prefer the same behavior. Figure |5] shows the case, where the 
individuals of different subpopulations dislike each other: 

(**(*)) : = ( J J ) , i-e-, {flU)) = [\ I ) • (46) 

This corresponds to a mutual influence of one subpopulation on the respective other. The 
computational results prove: 
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• The individuals avoid behaviors which can be found in the other subpopulation. 

• The subpopulation a = 1 with the lower indifference L\ < L 2 is distributed around 
the prefered behavior X\ and pushes away the other subpopulation! 

Despite the fact that the initial behavioral distribution P a (x, 0) agrees in both subpopu- 
lations, there is nearly no overlapping of Pi(x, t) and P 2 {x,t) after some time. This is a 
typical example for polarization phenomena in the society. 

In figure ||], we assume that the individuals of subpopulation 2 like the individuals of 
subpopulation 1 and, therefore, do not react on the behaviors in subpopulation 1 with 
avoidance processes: 

M) = (!?). (am) - ( o I ) • («> 

As a consequence, P2(x,t) remains unchanged with time, whereas Pi(x,t) drifts away 
from the prefered behavior x\ due to avoidance processes. Surprisingly, the polarization 
effect is much smaller than in figure |5|! The distributions Pi(x,t) and P 2 (x,t) overlap 
considerably. This is, because the slope of P 2 (x, t) is smaller than in figure [| (and remains 
constant). As a consequence, the probability for an individual of subpopulation 1 to meet 
a disliked individual of subpopulation 2 with the same behavior x can hardly be decreased 
by a small behavioral change. One may conclude, that polarization effects (which often 
lead to an escalation) can be reduced, if individuals do not return dislike by dislike. 



5 Empirical determination of the model parameters 

For practical purposes one has, of course, to determine the model parameters from empiri- 
cal data. Therefore, let us assume to know empirically the distribution functions P%(x,ti), 
[the interaction rates ^(iz)] and the effective transition rates w^(x'\x; t{) (x' ^ x) for 
a couple of times ti G {t , . . . ,t L }. The corresponding effective social fields V^(x,ti) and 
diffusion coefficients Q e ai j(x,ti) are, then, easily obtained as 

X 

V a c (x, U) :=- J dx' ■ K e a (x', t{) (48a) 

with 

:= jMM-XiMWxiU), (48b) 
h 

and 

Q e aij (x,ti) := Jd n x'(x' i -x i )(x' j -x j )w^(x'\x;t l ). (49) 
h 

Much more difficult is the determination of the utility functions U*(x,ti), U£(x,t{), the 
distance functions D e a (x', x\ t t ), and the rates z^(*z), ^(*z) := Kbi^i) filial), Vab(ti) '■ = 
u ab(ti)fab(ti)- This can be done by numerical minimization of the error function 

P e a(x,tl)\ , (50) 



a=l Z=0 i,i'ffl 



w*(x'\x;ti) 



D a (x',x;ti) 



g a (x',x;ti) 
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for example with the method of steepest descent PHI - In (|50|), we have introduced the 
abbreviation 



^(x'^-U^xfy) 
[" ab\"i' I ' M b V,"" i "l J 1 " ab\ u i- ) b '■) "l J J ^ 

6=1 

(51) 

It turns out (cf. [[HJ), that the rates f a (tz) have to be taken constant during the 
minimization process (e.g., u a {ti) = 1), whereas the parameters U a (x,t{), U a (x,ti), 
v lb(ti) '■= Kb^fabiti) and Vabfa) := v c ab {ti)fl b {ti) are to be varied. For l/D a (x',x;ti) 
one inserts 

1 n a (x',x;ti) 



D a (x',x;ti) d a (x',x;ti) 



(52a) 



with 



n a (x', x; U) := w a c (x'\x\ U)g a (x' ', x; U) [P a °(x, t/)J + w"(x\x'; U)g a (x, x'\ U) [i^(aj', tj) 

(52b) 

and 



d a (x',x;ti) := |3- a (a/,£c;^)P*(cc,^)J + [9a(x, x'; ti)P°(x', U) 
follows from the minimum condition for D a (x', x\ t{) (cf. ||19|| ). 



(52c) 



Since F may have a couple of minima due to its nonlinearity, suitable start parameters 
have to be taken. Especially, the numerically determined rates z^(tz) and ^(tz) have to 
be non-negative. 

If F is minimal for the parameters U a (x,ti), U a (x,ti), D a (x',x;ti), v a (ti), ^abi^i) and 
u ab(ti)i ^is is (as can easily be checked) also true for the scaled parameters 

U^x,^) := U a (x,tt) - C a (tt) , 

u:(x,ti) : = u'fatfi-crfa), 

D a (x',x;ti) 



Dl{x',x-ti) 
<{ti) 



< c b(ti) ■- 



In order to obtain unique results we put 



D a (U) 

Vg(tl) 

D a (t t ) ' 



(53) 



xefi 



xen 



and 



E 



'en x ] t[ 



E i, 



(54) 



(55) 
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which leads to 



C«(f,) := , C a (U) := , (56) 

xSQ x£fl 



and 



E 



DJx\ x\ ti 



Daiti) ■ E 1 



C a (ti) and C a (ti) are mean utilities, whereas D a (t{) is a kind of imzi o/ distance. 



(57) 



The distances D°(x',x;t) are suitable quantities for multidimensional scaling |42| , |43| . 
They reflect the "psychical structure" (psychical topology) of individuals of subpopula- 
tion a, since they determine which behaviors are more or less related (compatible). By 
the dependence on a, Dl(x',x;t) distinguishes different psychical structures resulting in 
different types a of behavior and, therefore, different "characters" . 



6 Summary and outlook 

In this article, a behavioral model has been proposed that incorporates in a consistent way 
many models of social theory: the diffusion models, the multinomial logit model, Lewin's 
field theory, the logistic equation, the gravity model, the Weidlich-Haag model, and the 
game dynamical equations. This very general model opens new perspectives concerning 
a theoretical description and understanding of behavioral changes, since it is formulated 
fully mathematically. It takes into account spontaneous (or externally induced) behavioral 
changes and behavioral changes due to pair interactions. Two important kinds of pair 
interactions have been distinguished: imitative processes and avoidance processes. The 
model turns out to be suitable for computational simulations, but it can also be applied 
to concrete empirical data. 



6.1 Memory effects 

The formulation of the model in the previous sections has neglected memory effects that 
may also influence behavioral changes. However, memory effects can be easily included 
by generalizing the BoLTZMANN-like equations to 

d r 

P a (x,t) = J dt' E [wlAx\x';t')P a (x',t') - w a t _ v {x'\x-t')P a {x^)\ (58a) 



with the effective transition rates 



w 



UiA^-t') :=™r'(aW)+EE E^V^'k^Onte,*'), (58b) 

b=i t/eQi/'en 



16 



and generalizing the Boltzmann-Fokker-Planck equations to 

t 



to 



n a 



d d 



2 /-^ dxi dxj 



Qtif(x,t')P a (x,t') 



with the effective drift coefficients 



d n x' (x't-xjw^ix'lx;?), 



(59a) 



(59b) 



the effective diffusion coefficients 



(59c) 



and 



wU'{x'\x;t') :=w t - t '(x'\x;l/) + J2 / d n y J d n y' wtf \x' ,y'\x,y;t')P b (y,t') . (59d) 

6=1 n a 



Obviously, in these formulas there only appears an additional integration over past times 
i! fill . The influence of the past results in a dependence of wf_ t ,(x'\x;t'), K 1 ^ 1 (x,t'), 
and Q l ~f (x,t') on (t - t'). The BoLTZMANN-like equations (f|) resp. the Boltzmann- 
Fokker-Planck equations (^) used in the previous sections result from (^) resp. (|59|) 
in the MARKOVian limit 



w%_ t ,(x'\x; t') := w a {x'\x; t)5{t - t') 
of short memory (where 5(.) is the Dirac delta function). 



(60) 



6.2 Analogies with chemical reactions 

The BoLTZMANN-like equations @ can also be used for the description of chemical re- 
actions, where the states x denote the different sorts of molecules (or atoms), and a 
distinguishes different isotopes or conformeres. Imitative and avoidance processes corre- 
spond in chemistry to self-activatory and self-inhibitory reactions. Although the concrete 
transition rates will be different from (|T0|) , (|lla|) in detail, there may be found analogous 
results for chemical reactions. Note, that the Arrhenius formula for the rate of chemical 



reactions jPJ can be put into a form similar to (|lla| ) |T9 
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Figure 1: Imitative processes in the case of one-sided sympathy and low indifference 
L a with respect to behavioral changes. 




Figure 2: As figure 1, but for high indifference L a with respect to behavioral changes. 




Figure 3: Imitative processes for mutual sympathy and low indifference L a in both 
subpopulations. 




Figure 4: As figure 3, but for high indifference L 2 in subpopulation 2. 




Figure 5: Avoidance processes for mutual dislike of both subpopulations. 




Figure 6: Avoidance processes for one-sided dislike. 



